Static properties of positive ions in atomic Bose-Einstein condensates 
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The excess number of atoms around an ion immersed in a Bose-Einstein condensate is determined 
as a function of the condensate density far from the ion. We use thermodynamic arguments to 
demonstrate that in the limit of low densities the excess number of atoms is proportional to the 
ratio of the atom-ion and atom-atom scattering lengths. For denser systems we calculate the excess 
number from solutions of the Gross-Pitaevskii equation using a model potential that has a 
attraction coming from the polarization of the neutral atoms and a hard core repulsion at short 
distances. We show that there exist in general many solutions to the Gross-Pitaevskii equation for a 
given condensate density, the maximum number of solutions being related to the number of bound 
states of the Schrodinger equation for the same potential. With increasing density, pairs of these 
solutions merge and disappear, implying a discontinuous change of the state of the system. 

PACS numbers: 03.75.Hh, 67.90. +z 



I. INTRODUCTION 

For the helium liquids, measurements on ions have 
served as a valuable probe of liquid properties. As exam- 
ples we may mention the use of ions to put in evidence 
vortex lattices in rotating liquid ^He and the measure- 
ment of ionic mobilities in superfluid ^He to elucidate 
scattering processes. The first experiments on ions in an 
ultracold gas of ^^Rb atoms were reported by the group 
in Pisa who produced ions by irradiating a rubidium 
condensate with laser pulses which ionize atoms through 
one- and two-photon absorption processes. 

Theoretically, the capture of atoms into weakly bound 
states of the atom-ion potential has been considered in 
Ref. • In this paper we consider the structure of a pos- 
itive ion in a Bose-Einstein condensate when there is no 
capture of atoms into bound states, and in particular we 
calculate the excess number of atoms associated with an 
ion. We shall demonstrate that this number is typically 
of order 10^, either positive or negative. 

The interaction between an atom and a positively 
charged alkali-metal ion (charge e), which are separated 
by a distance r, is given at large distances by the polar- 
ization potential caused by the electrostatic field £cs due 
to the ion, £cs — e/47reor^, which gives rise to a change 
in the energy of the neutral atom given by = 
where a is the polarizability of the atom. Writing the 
polarizability as a = Aneoa, where a has the dimension 
of volume, the energy shift of an atom due to the ion 
becomes 



(1) 



where Cq = e^/47reo. At short distances (r^lOag) the 
potential has a repulsive core. An important character- 
istic length, which we denote by /?4, may be identified 
by equating the kinetic energy /2m0l to the potential 



energy V{[3/C), resulting in 



„ , a m 

P4 = \ . 2 

ao nie 

Here m denotes the mass of a neutral atom, rric is the 
electron mass and oq = /nice^ ~ 0.53 A is the Bohr 
radius. Using the measured values a = 320 Aq for ^^Rb 



and a — 163 Oq fo 



r 23 Na, one finds I3f° « 7150ao and 
(3^^ K, 2620ao. The quantity (3^ gives the distance from 
the ion beyond which the zero-energy atom-ion wave 
function ceases to oscillate, and it sets the scale of atom- 
ion scattering lengths, but their actual values depend on 
the details of the potential at short distances. 

We begin by deriving from thermodynamics a general 
expression for the excess number of atoms around an ion 
and show that in dilute systems the excess number de- 
pends only on the ratio of the atom-ion and atom-atom 
scattering lengths. As we shall see, this approach sug- 
gests that the number of atoms associated with an ion is 
typically of the order of 10-100, but that it may be either 
positive or negative. In denser systems the excess number 
must be obtained from microscopic considerations, and 
we shall determine the structure of an ion immersed in 
a Bose-Einstein condensate at zero temperature, assum- 
ing that atom-atom interactions may be described within 
the framework of the Gross-Pitaevskii (GP) mean-field 
approach. We present solutions of the GP equation for 
a number of potentials which include a hard core repul- 
sion, an attractive square well, and one which resembles 
the atom-ion interaction, a hard core with a l/r'* attrac- 
tion at larger distances. 

For a given inner boundary condition, the Schrodinger 
equation has only one solution for a given value of the 
energy. By contrast, the GP equation, because it is non- 
linear, can have more than one solution for a given chem- 
ical potential. For potentials like the atom-ion one that 
can support two-body bound states, we shall find that 
at low densities there are 2vs + 1 solutions of the GP 
equation, where is the number of bound states of the 
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Schrodinger equation for the same potential. With in- 
creasing density, pairs of solutions merge and disappear 
until there is only a single solution with no nodes. We 
shall illustrate this behavior for two potentials, an attrac- 
tive square well and one with an attractive \jr'^ tail. An 
important question is which of these solutions is physi- 
cally relevant. At low condensate densities, one expects 
the wave function close to the ion to resemble the zero- 
energy solution of the Schrodinger equation, and to have 
I's nodes. This will be the case unless inelastic processes 
can populate lower-lying states. We find that with in- 
creasing condensate density, this solution ceases to exist. 
This indicates that the evolution of the state with den- 
sity cannot be continuous even in the absence of inelastic 
processes. 

The plan of the paper is as follows. In Sec^jwe present 
thermodynamic considerations. Section IIIII contains a 
description of the asymptotic behavior of the condensate 
wave function far from the ion. In Sec.|^we consider two 
simple potentials to illustrate important general features 
of our results, and in Sec. we analyse the case of a 
potential that, Hke the actual atom-ion potential, behaves 
as at large distances. We calculate the excess number 
of atoms from numerical solutions of the Gross-Pitaevskii 
equation for a given background condensate density. The 
concluding section. Sec. IVIL discusses our main results. 
In an appendix, we address the question of the validity 
of the Gross-Pitaevskii equation in the present context. 



II. THERMODYNAMIC CONSIDERATIONS 

We wish to calculate the excess number of particles as- 
sociated with an ion. To define this quantity precisely, 
we imagine adding an ion to a condensate. This will gen- 
erally change the density of atoms far from the ion by 
an amount that varies as 1/1^, where V is the volume of 
the system. A natural definition of the excess number 
of atoms AiV associated with the ion is the number of 
particles that must be added to keep the atom chemical 
potential constant, since this will ensure that the prop- 
erties of the condensate far from the ion are unaltered 
by the addition of the ion. In terms of the microscopic 
density of atoms n(r) around the ion, the excess number 
is given by 



AiV = 47r 



dr \n{r) — no] 



(3) 



where no is the density of atoms at large distances from 
the ion. 

This is analogous to what has been done earlier to 
calculate the excess number of ""^He atoms associated with 
a "^He impurity in Hquid ''He 0|. We shall denote the 
energy per unit volume as f (na,ni), where rta and n\ are 
the number densities of atoms and ions, respectively. The 
chemical potential of the atoms is given by 



and therefore the condition that this be unchanged by 
adding one ion and AiV atoms is 



AN = 0, 



AiV: 



(5) 



(6) 



When the density of ions is sufficiently low, d£/dni is 
equal to the energy change Ci when one ion is added to 
the condensate, and therefore 



AiV 



(7) 



One may also calculate AiV from the change AF in the 
thermodynamic potential F = E — fi^N when a single 
ion is added to the system at constant fia,. Here E is the 
total energy and N the total number of atoms. Since the 
number of atoms is given by 



N : 



it follows immediately that 



AiV = 



dF 
dAF 



(8) 



(9) 



Provided the volume considered is large compared with 
the scale of the atom excess around the ion, AF will be 
independent of the volume. 

Let us begin by making estimates for a dilute gas. Pro- 
vided the scattering of atoms by atoms and of atoms by 
ions may be treated as independent binary events, the 
energy density may be expressed in terms of the scatter- 
ing lengths associated with the atom-atom and atom-ion 
interactions. If ion-ion interactions are neglected, we may 
write 



f (na,ni) = ^C/aa"a + UaiU^ni, 

and therefore from Eq. © we obtain 



AiV = 



C/aa' 



(10) 



(11) 



The mean-field interaction constant Uji for species j and 
I, which may be either atoms (a) or ions (i), is related to 
the scattering length Uji by 



(12) 



where viji — riijmi/ [mj + mi) is the reduced mass of the 
two particles. Our result can therefore be expressed as 



AiV 



(13) 
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If, as in Ref. the ion is obtained by photoionization 
of the condensate itself, the latter expression reduces to 

AiV = -aai/flaa- (14) 

To obtain an order of magnitude estimate of the excess 
number of atoms associated with an ion, we note that the 
characteristic scale for the magnitudes of atom-ion scat- 
tering lengths loail is set by /34, given in Eq. l(2l, while the 
scale for the magnitudes of atom-atom scattering lengths 
Ittaal is set by 

1/4 



/36 



mo 



(15) 



Here Cg is the coefficient of in the van der Waals 
interaction, expressed in atomic units. Thus we arrive at 
the estimate 



IA7VI 



/ -2 \ 1/4 



(h 

/36 " \2m,Cfi 



(16) 



which is of order 35 for Rb and 25 for Na. 

The fact that the excess number of atoms is so large 
indicates that it may well be a poor approximation to 
regard the ion as a free particle, with mass equal to the 
bare ion mass. Rather, the recoil of the ion will be sup- 
pressed by the other atoms surrounding the ion, and if 
AiV 3> 1 it will be a better approximation to regard the 
ion as being stationary. In that case the excess number 
of atoms will be given by 

aai(TO) 



AiV : 



2aa 



(17) 



where the argument of Cai indicates that the scattering 
length is to be evaluated for a reduced mass equal to 
the atom mass. Expression IjlTII gives a value for AA^ 
that is typically of the same order of magnitude as that 
given by Eq. ifTfit . However, we stress the fact that the 
estimate for AA^ depends sensitively on the value of the 
effective mass of the ion, since the atom-ion potential 
has many bound states, and therefore relatively small 
changes in the reduced mass can result in large changes in 
the scattering length. Given that in the limit of low atom 
density the magnitude of the excess number of atoms is 
expected to be very much greater than unity, the result 
itTljl will generally not give a realistic estimate even in 
that case. 

The perturbation induced by the ionic potential is very 
strong. Therefore the question arises of whether the cus- 
tomary assumption of an essentially zero range for the 
atom-atom interaction is valid. We address this point 
in Appendix A, where we argue that the corrections to 
the GP result should not be large for the properties of 
interest here. 



tion in the vicinity of an ion extends to large distances 
from the ion and involves many atoms, we expect that 
the effective mass of an ion and its dressing cloud will 
be much larger than that of an atom, and we may re- 
gard the ion as being static. To describe the structure of 
the condensate in the vicinity of an ion we must therefore 
calculate the structure of the condensate in a static exter- 
nal potential given by the atom-ion interaction. Provided 
the length scale on which the condensate wave function 
V' varies in space is sufficiently large, we may do this 
by employing the Gross-Pitaevskii equation with the in- 
teraction of atoms with the ion included as an external 
potential. 



Co 1*1' 



(18) 



Here and in what follows we shall denote the chemical 
potential of an atom by /i, and for simplicity we have 
written Uq = i7aa — 47r?i^aaa/?n, since maa — to/2. We 
wish to find solutions that tend to a constant at large dis- 
tances from the ion, and since the potential is spherically 
symmetric, these solutions depend only on the radial co- 
ordinate r. Thus Eq. II18|I becomes 



2to dr"^ 



X = MX, 



(19) 



where x = r-f/;. 

The behavior of the condensate wave function at large 
distances depends on the nature of the potential V{r). 
On linearizing the GP equation lfT8|l and making use of 
the fact that the chemical potential is related to the con- 
densate wave function -i/iQ at large distances by the re- 
lation /i = noUo where no = IV'oP, one finds that the 
deviation 



(5-0 = i/) - i/'o 



(20) 



of the condensate wave function from its asymptotic value 
satisfies the linearized GP equation 



-—^^ + V{r) + 2Uon^]5ij^-V{r)^o. (21) 



For potentials with a finite range, one may neglect the 
potential at large distances from the ion, and the devia- 
tion that vanishes for r oo is thus given by 



(5-0 oc 



(22) 



III. MICROSCOPIC THEORY 

We now turn to microscopic considerations. Since, as 
we shall see, the distortion of the condensate wave func- 



where fc^ = \/2/ i and ^ is the healing length for the bulk 
condensate, 



1 



VSTTOaano 



(23) 
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For a potential, such as the atom-ion potential, that 
falls off at large distances less rapidly than the solution 
l|22ll . the behavior is different. The leading term in the 
solution for large r is then the Thomas-Fermi result V'tf, 
given by 



V^(r)+[/o|VTF|'=/i, 



(24) 



which, for the atom-ion potential with asymptotic form 
given by Eq. 1^, amounts to 



riTFir) = no 



Vir) 



Uo 



«o 1 + 



(25) 



or, to first order in V, 



'TF 



Vjr) 



(26) 



where we have taken V'o to be real. The density pertur- 
bation at large distances is seen to be always positive. 
Corrections to this result for smaller r may be calcu- 
lated from Eq. (I21II by neglecting the potential on the 
left hand side of the equation. The resulting differential 
equation may be written in terms of a function Sx defined 
by Sx — rStp, 



2m 

I^X^~j^rV{r)i^o- (27) 



This differential equation may be solved exactly in terms 
of exponential integrals, the two linearly independent so- 
lutions of the homogeneous equation being exp(— fc^r) 
and exp(A:^r). By inspection of lf27jl it is evident that 
the leading term for large r of the particular solu- 
tion to the inhomogeneous equation is given by Sx = 
— (2m / fi^ k'^)rV {r)ipo , which yields the Thomas-Fermi ex- 
pression Il2(i|l . The correction to this result may be ob- 
tained from the exact solution, but it is simpler to iterate 
Eq. lf27jl by moving the term d'^Sx/dr^ to the right hand 
side and replacing x in it by the Thomas-Fermi solution. 
This results in 



2m 



V{r) 1 + 



12 



(28) 



The leading correction to the Thomas-Fermi result for 
Sip given in Il2(i|l is thus seen to be proportional to r~^. 
Since we have already neglected the potential energy on 
the left hand side of 1I2TI1 . we cannot by this method 
obtain higher-order corrections to the particular solution 
than the one exhibited in H28II . 

By keeping in the general solution only the exponen- 
tially decaying term we thus get the asymptotic result 



Vjr) 
2noUo 



-kcr 



c- 



(29) 



where C is an arbitrary constant. 



For r ^ 00 the asymptotic behavior of the solution is 
always given by the TF result. However, whether or not 
this behavior is relevant for determining the structure of 
most of the cloud of atoms surrounding the ion depends 
on the relative size of the two characteristic lengths, /S^ 
and On the one hand, for P4 £, most of the cloud 
will be described by the TF approxmation, and only at 
distances less than ^ will the exponential term become 
important. On the other hand, for ^ ^ /34 (i.e. for low 
external density) the structure will be dominated by the 
exponential term, and the TF tail will become impor- 
tant quantitatively only at very large r. At shorter dis- 
tances from the ion, the mean-field energy becomes small 
compared with the atom-ion potential and the GP equa- 
tion reduces to a good approximation to the Schrodinger 
equation. 



IV. SIMPLE MODEL POTENTIALS 

Before presenting results for the attractive poten- 
tial we begin by examining two simpler model potentials, 
a repulsive hard-core and a spherical well. 



A. Hard-core potential 

Consider an interacting Bose-Einstein condensed gas in 
the presence of a repulsive hard-core potential of radius 
R. This model may be treated analytically in both the 
small and large core radius limits. The solution to the 
GP equation at large distances from the ion is given by 
Eq. JSl, 



1 + C 



exp(— fc^r) 



(30) 



If one assumes that this expression holds for all r greater 
than R, we can determine the constant of proportionality 
C by imposing the boundary condition ipiR) — 0. This 
gives C — —Re^^^y/n^. For r close to R this has the form 
of the scattering solution for the Schrodinger equation, 
ip ^ 1 — R/r. In fact, for i?/^ <^ 1 this solution becomes 
essentially exact, since this function fails to satisfy the 
GP equation only in the region where r ~ i?, and in 
this region the total change in the slope dx/dr of the 
radial wave function is small and may be neglected. As 
an illustration of this fact, we calculate the excess number 
of particles, which is given by Eq. (pj, and find 



AiV = Anno 



dr 



-2R- 



31^ 
2V2 



-fce(r-i?,) 



i?3 

3 



R' 



-2fcj(r-fl) 



AttuoR^ 



(31) 



5 



For ^ ^ i?, this reduces to 



AiV = -■ 



_R_ 



(32) 



Let us now compare this result with the one derived on 
the basis of thermodynamic arguments. For a hard-core 
potential the scattering length coincides with the core 
radius. Since we have assumed the ion to be station- 
ary, its effective mass is taken to be infinitely large, and 
therefore the reduced mass for the ion and an atom is m, 
rather than the value m/2 one obtains for an ion and an 
atom with equal masses. Thus, this result is in precise 
agreement with Eq. ifT^jl . 

When the core radius is much larger than the healing 
length, the wave function reaches its asymptotic value 
on a length scale that is short compared to R. In Eq. 
ifTflll we can therefore replace the factor appearing in 
the nonlinear term by the constant and we are left 

with an effectively one-dimensional GP equation whose 
solution is 



r - R 

4' = \Aio tanh — -= — , r > R, 
v2^ 



(33) 



and zero otherwise, as may be seen by inspection. The 
excess number of particles is given by 



4 i?^ 
AN = --nR^no 



\/2aaa? 



(34) 



where the leading term is due to exclusion of atoms from 
the core. 



B. Attractive square well 

We next consider a more physical potential, an attrac- 
tive well: 



y(r) 



2m ' 



r < R, 



(35) 



V{r) — otherwise. Like the actual ion-atom potential, 
this can have bound states for the two-body problem. 
With this potential we shall be able to examine how so- 
lutions of the GP equation disappear as the condensate 
density increases. The GP equation H19|l reads 



47raaa I ^ no I x = 

(36) 

(where x = rip and 9{x) is the step function) and the 
scattering length for this potential is 



R 1 



tan k^R 
koR 



(37) 



Since this equation is nonlinear, there can be multi- 
ple solutions for the same boundary conditions (i.e. the 



same bulk density no). As we will show in the follow- 
ing, for small no it has 2i^s + 1 solutions, where vs is 
the number of nodes of the zero-energy Schrodinger so- 
lution ips or, equivalently, the number of bound states of 
the Schrodinger equation. In the low background den- 
sity limit, inside the well the solution with the maximum 
number of nodes approaches ips, i-e. ip{r) oc sm{kor)/r, 
while outside it tends towards the uniform density no 
with the asymptotic behavior given in Eq. II22|I . Sip{r) cx 
exp{—k^r)/r. 

With increasing no, the mean-field repulsion between 
the atoms makes the effective potential shallower, which 
tends to push nodes of the wave function outwards. At 
the same time, the increase in the chemical potential has 
the opposite effect on the nodes. What we find is that if 
the zero-energy solution of the Schrodinger equation has 
vs nodes, for low condensate densities the GP equation 
has one solution with no nodes, and two solutions with 
any nonzero number of nodes less than or equal to ■ 

To demonstrate this, we analyze separately the be- 
havior of the wave function inside and outside the well, 
and match them at some intermediate point, which for 
this particular potential we take to be the edge of the 
well. Specifically, we integrate out from the origin, where 
X = 0, for different choices of the derivative of x at r = 
and calculate Tp and Tp' at the boundary r = R. These 
trace a curve m ip — ip' space. Then we integrate in- 
wards from large distances, where the solution is defined 
by the proportionality constant C of the Yukawa asymp- 
totic form, Eq. \22l . As C is varied, another curve in 
ip — tp' space is traced out. If the mean-field interac- 
tion could be neglected for r < R, the ratio ip'{R)/'p{R) 
would not depend on the normalization of the wavefunc- 
tion, and therefore the curve corresponding to the inner 
boundary would be a straight line through the origin. In 
the presence of atom- atom interactions, the ratio tp' /tp 
obtained by integrating outwards traces out a spiral. For 
low no this crosses the ip' axis a number of times equal 
to the number of nodes the zero-energy solution of the 
Schrodinger equation has inside the potential. This fol- 
lows from the observation that for low x'(0) the solution 
will have the same number of nodes inside the potential 
as the zero-energy solution of the Schrodinger equation, 
while for very large values of x'(0) the effects of the mean 
field will be so strong that the solution has no nodes in- 
side the potential. 

The corresponding plot obtained by integrating in- 
wards has two branches, depending on whether ip{r ^ 
oo) is positive or negative. Examples of the plots are 
given in Fig. for parameters such that i^s = 3. For 
low no, there are 2i/s -I- 1 intersections of the two sets of 
curves, corresponding to solutions of the GP equation. 
This is illustrated in Fig. Q^. As no increases, pairs of 
solutions with the same number of nodes merge and dis- 
appear, as shown in Fig. ^p. Eventually, at sufficiently 
high values of no only the nodeless solution survives. 

In Figs, m and 121 we show how, increasing the external 
density, the solutions with the highest number of nodes 
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Figure 1: (Color online) Behavior of ipiR) and ip'iR) for the 
solution inside the well (solid line) and outside it (dashed and 
dot-dashed lines for ip{r oo) = ii/'o, respectively). In the 
plots we have set fco^ = 9, which gives three bound states 
for the Schrodinger equation. We measure energies in units 
of /l^fco/Sm and lengths in units of R. The calculations were 
performed for Uo = 0.45 in these units, but results for other 
values of Uo may be obtained by scaling, since for a given 
chemical potential, ij) and tp' vary as Uq^^"^ . The symbols near 
intersections indicate the number of nodes of the solution, and 
for this case vs — 3. The upper panel (a) is for = 0.45, and 
the lower one (b) for /i = 2.9, just above the value /i = 2.52 
at which the two solutions with 3 nodes merge and disappear. 



actually merge. For densities higher than this critical 
value, the only solutions are ones with a smaller number 
of nodes. 

Despite its short-range character, the model given 
above captures the main features of the solutions of the 
GP equation for the atom-ion potential, which is long- 
ranged. We note that the discontinuous behavior does 
not occur in a one-dimensional model. 




1 2 3 4 5 

r/R 




1 2 3 4 5 

r/R 



Figure 2: (Color online) Two solutions of the Gross-Pitaevskii 
equation for the attractive square well potential. koR and Uo 
are the same as in Fig. Q For the upper panel, the chemical 
potential is 0.45, as in Fig. [^i, while for the lower one it is 
2.5, just below the value at which the solutions merge. The 
solutions both have three nodes, and are the first to merge as 
the chemical potential increases. 



V. THE r" POTENTIAL 

We now turn to a more realistic potential with the 
same r""' behavior as the actual atom-ion interaction 
at large distances. For definiteness, we consider param- 
eters appropriate for a ®^Rb condensate, and we take 
Qaa = 100 flQ. At large distances, we take the atom-ion 
potential to be given by Eq. Q with a = 320aQ. The 
wave functions are sensitive to the short-range behavior 
of the potential, but we may obtain illustrative results 
by cutting the potential off by a hard core of radius 
R. Since many atoms are bound to the ion, we assume 
the ion to be static and set m^i = to. The atom-ion scat- 
tering length of such potential may be calculated in the 
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WKB approximation, and is given by 

flai = Pi cot 



(38) 



The number of bound states allowed by the potential can 
be estimated by increasing the potential strength from 
to its actual value. A bound state appears each time the 
scattering length diverges, and therefore the number of 
bound states is given by 



ttR 



(39) 



where Int(a;) denotes the integer part of x. To model 
actual atom-ion potentials, a physically reasonable value 
of -R would be ~ IOoq. However, the properties of the 
wave function of most importance here are those at rel- 
atively large distances, r^lO'^oo, so we take R — SOOoo, 
since this should give us the correct physical behavior for 
the distances of interest. We do not expect the qualita- 
tive behavior of the wave function to depend on R, even 
though quantities Hke the scattering length do, and we 
have verified this numerically. 

We now describe numerical solutions of the GP equa- 
tion that approach a constant density uq far from the 
ion. Just as for the finite-range potential considered in 
the previous section, there is generally more than one so- 
lution for a given value of the chemical potential, and for 
small external densities one expects 2z^s + 1. In Fig. we 
show the wave functions corresponding to the two states 
with the highest number of nodes, namely seven for the 
parameters chosen, in agreement with the quasiclassical 
result Ij39|l . The free energy, for a given condensate den- 
sity no, is highest for the states with the highest number 
of nodes, and decreases as the number of nodes decreases. 

In the absence of inelastic processes, we expect only 
the uppermost state of the ionic potential to play an im- 
portant role in the capture process, since it is the only 
one with an appreciable overlap with the continuum wave 
function representing the unbound atoms 

The excess number of atoms is given in terms of the 
atomic density distribution by Eq. (pj or, alternatively, 
from the free energy F = E—fiN by Eq. JH}. In Fig. 01we 
show results obtained from our numerical simulations by 
both methods. In the Hmit of very low condensate density 
we get values for A7V in accord with the thermodynamic 
arguments in Sec. |n| The consistency of the two methods 
of calculation has been confirmed for core radii that give 
scattering lengths in the range |aai| < 5000 oq. 

The figure shows that the excess numbers of atoms for 
two states with the same number of nodes become equal 
at the density above which the solutions no longer ex- 
ist. This is to be expected, since the solutions become 
identical at this point. In Appendix A we use quasiclassi- 
cal arguments to estimate the density at which solutions 
merge and disappear, and these are in good agreement 
with the numerical results. 

In the detailed calculations described so far we have 
focused attention on states with close to the maximum 
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Figure 3: (Color online) Condensate wave functions for the 
two uppermost states in the 1/r* potential with the param- 
eters given in the text for no = lO^'^cm^'^. Both states have 
seven nodes, but the resolution of the figure is inadequate to 
exhibit the rapid oscillations for r close to R. The state that, 
in the dilute limit, becomes the zero-energy solution of the 
Schrodinger equation is given by the solid line. 



number of nodes. In particular, in the low-density limit 
and in the absence of inelastic processes that can cause 
the system to relax, one would expect the state of the con- 
densate to be the one that close to the ion resembles the 
zero-energy solution of the Schrodinger equation. How- 
ever, three-body processes can relax the system, thereby 
populating states with lower numbers of nodes. To cal- 
culate properties of such a system, one could start with 
a many-particle wave function of the Hartree-Fock type 
in which more than one single-particle state is occupied, 
and solve the Hartree-Fock equations. This is, however, 
beyond the scope of this paper because the density of 
atoms rises to values sufficiently high that the dilute gas 
approximation for the interaction energy employed in the 
GP approach fails at relatively large distances from the 
ion. To see this, we note that the density of atoms far 
from the ion will be given by the Thomas-Fermi approx- 
imation, Eq. Ij25ll . The dilute gas approximation is valid 
provided njaaap <C 1. This condition becomes 



V{r) 



87rr4 



< 1 



or 



r» (/34|aaa|)'/V2, 



(40) 



(41) 



which for rubidium (a^a ~ lOOao) implies that the GP 
equation is valid only for r ^ 400ao for such states. 



VI. CONCLUSIONS AND DISCUSSION 

In this paper we have investigated solutions of the 
Gross-Pitaevskii equation for a Bose-Einstein conden- 
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Figure 4: (Color online) Excess number of atoms around a 
single ion as a function of the bulk density. The dashed line is 
the dilute limit appropriate for a fixed ion, AA'^ — — aai/2aaa 
(R = 300 ao gives Oai ~ — 1980ao for an infinitely massive 
ion). Results are shown for the four uppermost levels for this 
potential (i.e. the two with 7 nodes and the two with 6 nodes, 
indicated respectively by the solid and dot-dashed lines). The 
lines are obtained from Eq. lO, the circles from Eq. lO. The 
inset exhibits the behavior at lower densities. 



sate in the presence of a positive ion. We find that for 
low condensate densities, there are 2z^s + 1 solutions for 
a given condensate density, where i^s is the number of 
bound states of the Schrodinger equation. With increas- 
ing condensate density, pairs of states become degenerate 
and disappear, and the state of the system must change 
discontinuously. An interesting challenge is to find ex- 
perimental evidence for such a behavior. 

We have calculated the excess number of atoms around 
an ion, and find that for the state that resembles the zero- 
energy solution of the Schrodinger equation it can be 
either positive or negative, depending on the sign of the 
atom-ion phase shift, and a typical magnitude is of order 
~ 10^ . The spatial size of the density disturbance around 
an ion is set \>y ji^ ^ \ fim. Our estimates indicate that 
the Gross-Pitaevskii equation should give a reliable first 
approximation for the wave function of such states. For 
states with fewer nodes, the density of atoms may reach 
values high enough that the GP equation fails. 

There are many outstanding problems. In most of the 
calculations we have assumed that the state of interest 
is that with the maximum possible number of nodes. 
More study is needed of inelasticprocesses that will cause 
atoms to relax to lower states Experimental studies 
will be valuable in providing guidance for future work. 
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Figure 5: (Color online) Difference in free energy for the 
states given in the previous figure: the solid lines are for the 
two states with 7 nodes and the dot-dashed line for one of the 
states with 6 nodes. The free energy is measured in units of 
W~^h^ /maQ. The other state with 6 nodes lies much lower, 
at around AF ^ —3 ■ 10~''fi^/mao. 
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Appendix: VALIDITY OF APPROXIMATIONS 

Here we examine the conditions under which it is a 
good approximation to replace the effective atom-atom 
interaction by the standard expression lfT2|l . Later we 
shall estimate the density below which the dilute gas re- 
sult l|i;^|l would be expected to hold. For nonzero wave 
numbers k, the quantity that enters the expression for 
the energy shift is —S/k, where 6 is the s-wave phase 
shift, rather than a Since a typical energy scale for 
changes in S/k is set by fi.^/2m/?|, while the potential 
depth is given by aeQ/2r'*, we expect that the scattering 
length approximation will fail when 

The Gross-Pitaevskii approach should therefore be valid 
if the phase shift due to the region where r <C f = 
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(/34/96)"'^^^ is negligible. To estimate this phase shift, we 
make a semiclassical approximation to the GP equation. 
This gives for the total accumulated phase out to a dis- 
tance r 



$(r) 



dr'\l'^[ti^Vir')^n{r')Uo]. 



(A.2) 



Deep in the ionic potential, the wave function is given to 
a good approximation by the semiclassical result, which 
has an amplitude 



l/2x 



(A.3) 



where Pdir) — [2mF(r)]^/^ is the classical momentum of 
a particle of zero total energy in the presence of the po- 
tential. For the r~* potential, the amplitude of the wave 
function is therefore independent of r, and therefore we 
may replace the mean-field energy to a first approxima- 
tion by a constant nUo, where h is independent of r. 
Expanding expression ljA.2|l in the deviation no — we 
find 



by an amount 



5$(f) « (no - n)aUPiPiy^^- 



(A.5) 



If we take the interior density to be of the same order of 
magnitude as that far from the ion, one finds 



(5$(f) 



(A.6) 



where the healing length is defined in Eq. H23|l . Since 
under experimental conditions the healing length is typ- 
ically comparable to P4, while is two orders of magni- 
tude smaller, this shows that the region close to the ion 
where the Gross-Pitaevskii equation fails is likely to be 
unimportant. 

On the basis of the above calculation, we may also 
estimate the density below which the low-density result 
11 3p is valid. Using the approximations above, we find 
that the total accumulated phase out to a distance ~ ^4, 
where the semiclassical treatment fails, is of order 



$(r) w <I>o(r) -I- (no - n) 



2mU( 



dr' - 



1 



/o 2y^-2mV{r')/h'^ 
(A.4) 

Due to the mean-field interaction, the accumulated phase 
out to a distance r ^ f = (/34/36)^^^ is therefore changed 



PI 



(A.7) 



This indicates that changes to the accumulated phase can 
be significant under typical experimental conditions. 
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